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ABSTRACT 

We have studied solar-like oscillations in ~ 800 red-giant stars using Kepler long-cadence photome- 
try. The sample includes stars ranging in evolution from the lower part of the red-giant branch to the 
Helium main sequence. We investigate the relation between the large frequency separation (Av) and 
the frequency of maximum power (i/ max ) and show that it is different for red giants than for main- 
sequence stars, which is consistent with evolutionary models and scaling relations. The distributions 
of i/ max and Av are in qualitative agreement with a simple stellar population model of the Kepler 
field, including the first evidence for a secondary clump population characterized by M > 2 M Q and 
y max ~ 40—110 //Hz. We measured the small frequency separations Sv 02 and Sv i in over 400 stars and 
duos in over 40. We present C-D diagrams for I = 1, 2 and 3 and show that the frequency separation 
ratios Sv 02 / Av and 8v m / Av have opposite trends as a function of Av. The data show a narrowing of 
the I — 1 ridge towards lower v mayil in agreement with models predicting more efficient mode trapping 
in stars with higher luminosity. We investigate the offset e in the asymptotic relation and find a clear 
correlation with Ai/, demonstrating that it is related to fundamental stellar parameters. Finally, we 
present the first amplitude-f max relation for Kepler red giants. We observe a lack of low-amplitude 
stars for i/ majt > 110 /xHz and find that, for a given i/ max between 40 — 110 //Hz, stars with lower Av 
(and consequently higher mass) tend to show lower amplitudes than stars with higher Av. 

Subject headings: stars: oscillations — stars: late-type 



1. INTRODUCTION 

Stars with convective envelopes show solar-like oscil- 
lations that are sensitive to th e physical processes gov- 
erning their interiors (see, e. g . , IBrown fc Gilliland! Il99i 
IChristensen-Dalsgaardl [2004D . Following the success of 
helioseismology, the detection of such oscillations in a 
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variety of stars holds great promise for improving our 
understanding of stellar structure and evolution. 

The traditional goal of asteroseismology is the accu- 
rate measurement of individual mode frequencies, which 
can be used to test stellar physics by comparing them to 
frequencies predicted by models. In a more general ap- 
proach, global oscillation parameters can be used. These 
include the average frequency separations, which are di- 
rectly related to properties of the sound speed in the stel- 
lar interior and therefore to fundamental stellar parame- 
ters. Other parameters are the amplitude and central fre- 
quency of the oscillation envelope, which are important 
for understanding the physics of driving and damping 
of these modes. The measurement of these parameters 
presents a valuable addition to classical methods such as 
spectroscopy, and a powerful tool to systematically study 
stellar evolution when oscillations are detected in a large 
ensemble of stars. 

Compared to main-sequence stars, red giants pulsate 
with larger amplitudes and longer periods, therefore re- 
quiring less sensitivity but longer and preferably con- 
tinuous time series for unambiguous detections. The 
first attempts to detect oscillations in G and K giants 
were focused on nearby targets such as Arct urus (jSmithl 
[19831 iCochranl [l988t iBelmonte etldl [l990l) . This was 
followed by campaigns targeting single stars and clus- 
ters using precise gr o und-based Doppler sp ectroscopy 
dFrandsen etafl 120021: iDe Ridder etafl l2006t ) and pho- 
tometry (|Stello et al I l2007r i. as wel l as using space-based 
photometers such as the HST ([Edmonds fc Gilliland! 
[1991: IGillilandl [20081 IStello fc Gillilandl 12009ft . WIRE 
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(iBuzasi et al.l2000HRetter et al.l2003tlStello et al]l2008l) . 
MOST (iBarban et al|200rtlKariinger et alJ^OOSa Mand 
SMEI (jTarrant et alj|20Q7ji . 

A breakthrough was achieved by CoRoT with unam- 
biguous detections of radial and non-radial modes in 
~ 800 red giants (iDe Ridder et al.1 [200l IHekker et all 
20091 ICarrier et al.ll2010D ~ With a maximum time series 
length of 150 d but limitations due to the low-Earth orbit 
of the satellite, the CoRoT detections were focused on 
low- mass He-core burning sta rs (the red clump). Using 
these data, iMiglio et al.l (|2009f ) performed the first popu- 
lation study using global oscillation parameters and con- 
cluded that the distributions were qualitatively in agree- 
ment with the c urrent picture of the star formation rate 
in our galaxy. iMosser et al.l (|2010l ) subsequently used 
a larger number of detections in the CoRoT sample to 
investigate correlations of various oscillation parameters. 

A new era of "ensemble asteroseismology" was re- 
cently entered with the laun ch of the Kepler space tele- 
scope (jGilliland et al.l l2l)H)i. First results for red gi- 
ants were based on 34 d of data (Bedding et al. 2010b; 
see also iStello et al.l l2010t IHekker et al.l l2010bl ) which 
demo nstrated the enormous potential of Kepler data. 
While Bedding ct al. (2010b) focused on a sample of low- 
luminosity giants, the goal of this paper is to systemati- 
cally investigate global oscillation parameters in the com- 
plete Kepler red-giant sample using data spanning up to 
138 d. We refer to our companion papers for the compar- 
ison of global oscillation paramete rs derived using differ- 
ent methods (jHekker et al.ll2"010d ) and t he asteroseismic 
determ ination of stellar masses and radii ([Kallingc r et al.l 
l2010bTl . 

2. OBSERVATIONS AND DATA ANALYSIS 

The Kepler space telescope was launched in March 
2009 with the principal science goal of detecting Earth- 
like planets around solar-like stars through the observa- 
tion of photometric transits. Kepler employs two ob- 
servation modes, sampling data either in 1 min (short- 
cadence) or 29.4 min (long-cadence) intervals. For our 
study of puls ations in red giants we used Kepler long- 
cadence data (| Jenkins et al.ll2010l ). which have a Nyquist 
frequency of 283 ^iHz. 

Kepler is located in an Earth-trailing orbit with space- 
craft rolls performed at quarterly intervals to redirect so- 
lar panels towards the Sun. Data are consequently sub- 
divided into quarters, starting with the initial commis- 
sioning run (10 d, Q0), followed by a short first quarter 
(34 d, Ql) and a full second quarter (90 d, Q2). The basis 
of our study consists of 1531 light curves for which data 
from all these quarters are available (total time span of 
137.9 d). We did not include stars thought to be mem- 
bers of the open clusters in the Kepler field (NGC6791, 
NGC6811, NGC6819 and NGC6866). 

Before we extracted asteroseismic information from the 
time series, two main instrumental artifacts had to be 
addressed. Firstly, in some cases data from different 
quarters show intensity discontinuities that are mainly 
caused by pixel shift s after reorientation of the spacecraft 
(jJenkins et al.lf20To| ) . Additionally, two safe-mode events 
in Q2 caused intensity drifts due to thermal effects, af- 
fecting in total ~ 7.5 days of the data. We discarded 
data affected by the safe-mode events and corrected for 
intensity jumps by robust fitting and correcting second- 
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Fig. 1. — Top panel: Kepler light curve of the oscillating red 
giant KIC5515314. The different quarters of data used in this 
work are indicated. Middle panel: Amplitude spectrum of the light 
curve shown in the top panel. Note that the ordinate shows the 
square root of power density, where the latter is power multiplied 
by the effective length of the dataset (calculated as the inverse 
of the area under the spectral window). Bottom panel: Spectral 
window, shown at the same abscissa scale as the middle panel. 
Note that the height of the peak at zero frequency is 1. 



order polynomials to QO, Ql and five separate segments 
of Q2 which were unaf fected by intensity jumps (see also 
iKallinger et alJfeOlObl ). Finally, we performed a simple 
4-cr outlier clipping to remove remaining outlying data 
points (in most cases < 0.1% of the total number of data 
points). 

The top panel of Figure [T] shows a typical light curve 
after the above corrections have been performed. The 
different quarters contributing to the final light curve 
are indicated. The main gaps in the datasets span 1.5 d 
between QO and Ql, 4.5 d between Ql and Q2, and 5 d as 
well as 2.5 d during the aforementioned safe- mode events 
in Q2. The middle panel displays the amplitude spec- 
trum of the light curve, showing the clear signature of 
solar-like oscillations with regularly spaced peaks cen- 
tered around 130 /iRz. The spectral window shown in 
the bottom panel of Figure [1] demonstrates the nearly 
continuous sampling of Kepler data. The weak peaks 
at multiples of 3.96 fiHz are caused by the angular mo- 
mentum dumping cycle of the spacecraft causing one re- 
jected datapoint every 2.9 d. At a level of less than 1% in 
amplitude, these aliasing artifacts are negligible for the 
analysis considered in this paper. 

To extract oscillation parameters from the corrected 
time series, several methods have been employed 
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(Hekkcr c taLll2010dt IHuber et al.l [20091 iKallinger et al 



2010U iMathur et all 120101: iMosser fc Appourchaux 
2009). For a detailed comparison of the result s of thes e 
different methods we refer to iHekker et all (|2010cf ). 
Unless otherwise mentioned, all results shown here 
are based on the method by IHuber et al.l ((2009). All 
values from that method were cross-checked with the 
others and we only retained stars for which at least one 
other method yielded consistent results. Uncertainties 
reported in this paper are based on calibrations of 
extensive (~ 50000) simulations of artificial time series 
with identical sampling to the Kepler data, including 
mixed modes, background noise and varying mode 
lifetimes between 15-75 d. 

Before reporting our results, we provide a brief sum- 
mary of the global oscillation parameters used in this 
work, their physical interpretation and the principal 
methods with which they were determined. 

3. GLOBAL OSCILLATION PARAMETERS 

3.1. Frequency of maximum power (v max ) 

As first argued by iBrown et al.l ([1991D , u max for sun- 
like stars is expected to scale with the acoustic cut-off fre- 
quency and can therefore b e related to fundamental s tel- 
lar parameters, as follows (Kjcldscn & Bedding 
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As shown by IStello et~a l. (2009) for stellar models and 
as observed in stars with well-determi ned fundamental 
parameters (see, e.g.. IStello et al.ll2008D . this scaling re- 
lation also holds for red giants. In these stars, z/ roax is 
a good indicator of the evolutionary stage, and ranges 
from ~ 10 /iHz for high-luminosity giants to ~ 250 /iHz 
for H-shell-burning stars in the lower part of the red- 
giant branch. It is measured by determining the max- 
imum of the power excess after heavily smoothing over 
several orders or by fitting a Gaussian function to the 
excess power. 

To illustrate our sample, Figure [2] shows all stars 
in a plot of ^ max versus T c g. Effective temperatures 
have been taken fro m the Kepler Input Catalog (KIC) 
(jLatham et al.l [20051 ) and we ha ve overlaid solar-scaled 
AST EC evolutionary tracks ([Christensen-Dalsgaardl 
120081 ) . In addition to the red giants, we also show a sam- 
ple of main-sequence and sub-giant stars for which oscil- 
lations had been d etected prior to the Kepler mission (see 
IStello et al.l 120091 and references therein). Figure [5] can 
be viewed as an asteroseismic H-R diagram in which, in 
the absence of parallaxes, we have used 1 /V ma x instead of 
luminosity. As can be seen from the evolutionary tracks 
in Figure [2j our sample spans a total range of masses of 
approximately 1-3 Mq , with temperatures ranging from 
4200 to 5200 K. For a detailed study of the fundamen- 
tal parameters of the stars in the sample we refer to our 
companion paper (IKallinger et alJl2010bD . 

3.2. Frequency separations 

According to the asymptotic relation fo r modes of low 
angular degr ee I and high radial order n (|Tassoullll980l ; 
Gough 1986), frequencies of solar- like oscillations can be 
described by a series of characteristic separations. Obser- 
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Fig. 2. — fmax versus effective temperature for all red giants 
in our sample ( diamonds), as well as main-sequence and sub-giant 
stars studied by Stcllo et al. (2009) (squares). Error bars have been 
omitted for clarity. The grey lines show soiar-met allicity ASTEC 
evolutionary tracks ( Christcnscn-Dalsgaard 2008) for a range of 
masses (note that the 0.8 Mq track has been evolved beyond the 
age of the universe). The dashed line marks the approximate po- 
sition of the red edge of the instability strip. 



vationally, these separation s can be expressed as follows 
([Bedding fc KieldserilMoh : 



v n .i = Av(n + -I + e) - 5v m 



(2) 



Here, Av denotes the mean large frequency separation 
of modes with the same degree and consecutive order. 
Av is directly related to the sound travel time across 
the stellar di ameter and probes the mean stellar density 
(|Ulrichlll986D . This means that Av is expected to scale 
as follows: 



Av 



(A//M Q )»- 5 (T cff /T cff , Q ) 
(i/io) - 75 



-Av Q . 



(3) 



In Equation ([2]), Svqi denotes the small frequency sepa- 
rations of non-radial modes relative to radial modes, as 
follows: 



<^01 = 77(^,0 + ^n+l,o) — v n,\ , 



&>03 = \(ynfl + fn+l,o) 



(4) 
(5) 



(0) 



Following iBedding et all (|2010bD , we have used Svqz in- 
stead of the more commonly used 8vx^ due to the broad- 
ening o f the I = 1 ridge in red giants caused by mixed 
modes ([Dupret et al.l 120091: iDeheuvels et al.l I2010D . For 
main-sequence stars, small frequency separations are sen- 
sitive to variations of the sound speed gradient near the 
stellar core, which changes as the star evolves due to the 
increase in its mean molecular weight. 

The phase constant e in Equation ([2]) has contributions 
from the inner and outer turning point of the modes. 
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This is usually expressed as e = k + a, where a is the 
contribution from the outer turning point, which is deter- 
mine d by the properties of the near-surface region of t he 
star (iChristens cn-Dalsgaar cTfe Perez HernandezHl992t ). 

3.3. Maximum mode amplitude 

The maximum mode amplitude is related to the 
turbulent convection me chanisms that exci te and 
damp the oscillations. iKjeldsen fc Beddingl (1995) 
fou nd that model predictions of sola r -like oscillations 
by IChristensen-Dalsgaard fc Frandsenl (|1983l ) implied a 
scaling for velocity amplitudes of 



(7) 



with s — 1.0 (see also lHoudek et al.1119991 ) . More recent 
calcul ations for main-sequence stars by iSamadi et al.l 
( 2007' ) have favored an expon ent of s = 0.7. 

IKieldsen fc Bedding! M l also argued that the oscil- 
lation amplitude A\ observed in photometry at wave- 
length A is related to the velocity amplitude: 



Ax oc — T eS 



(8) 



wher e r = 1.5 if the osci l lation s are adiabatic. How- 
ever, [Kicldsen & B eddingl (|1995h found a better fit to 
observed amplitudes in classical pulsators with an expo- 
nent of r = 2.0. 

Observationally, iStello et al.l (|2010D found agreement 
with s = 0.7 (assuming r = 2.0) for cluster red giants 
observed with Kepler, while IMosser et al.l ()2010D found 
a best fitting value of s = 0.89 ± 0.02 (assuming r = 1.5) 
for CoRoT red giants. As for f max , mode amplitudes 
are usually determined by heavily smoothing the power 
spectrum to eliminate variations caused by the stochastic 
nature of the signal (jKieldsen et al.l 120081 ) or by fitting 
a Gaussian function to the power excess envelope. In 
both cases, it is important to determine an accurate fit 
to the stellar background contribution. To account for 
the effects of averaging from the 29.4 min integrations, 
we divided the observed amplitude with a sine function: 

A rea l = A obs /sinc(f 
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Fig. 3. — Upper panel: Au versus v max for the Kepler red giants. 
The fitted relation is indicated by the red dashed line, with red dot- 
ted lines marking 3- ff uncertainties. T he blue dashed-dotted line 
shows the relation b y[Stcllo ct al. (2009). Note that the relation by 
IMosser" et al.l pOfOh and the average relation for Kepler stars de- 
rived using seven different analysis methods are indistinguishable 
from the red dashed line. Lower panel: Residuals after subtracting 
the fit indicated by the red dashed line. 



coefficients in Equation (j9]), we performed least-squares 
fits to the results of seven different methods. Note that 
we have omitted all stars with f max > 230 /zHz from the 
analysis of the Av-v max relation due to the difficulty of 
determining accurate ^ max values close to the Nyquist 
frequency. As can be seen in Table 1, there is good agree- 
ment between the r esults of different me thods, as well as 
with the result by IMosser et al.l (|2010f ) that was based 

on CoRoT stars with v m ax < 100 fjR z. 

As already shown by IStello et ail (|2009l ). a value of 
/3 ~ 0.75 is expected from the scaling relations for u max 
and Av. Raising Equation ([T]) to the power of 0.75 and 
dividing by Equation ([3]) yields 



4. RESULTS 

4.1. The Av - ^ max relation 

The growing number of detections of solar-like oscilla- 
tions across the HR diagram has revealed a tight power- 
law relation between Av and ^ max : 



Av = a(^ max /^Hz) 
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(9) 



Using a fit to r ed giants and main-sequence stars, 
IStello et all ([2009h derived a = 0.263 ± 0.009 //Hz and 
P = 0.772 ± 0.005, noting that models predicted slightly 
different relations for red giants than for main-sequence 
stars (see their Figure 3). The relation wa s also seen for 
the C oRoT red giants (i/ max < 100 mHz) by IHekker et abl 
(2009) and subsequently by IMosser et al.l (<2010D . with 
the latter finding a = 0.28±0.02 /iRz and /3 = 0.75±0.01. 

The Kepler data allow us to investigate this relation 
for the first time over a large part of the red giant branch. 
Figure |3] shows the relation for our sample. To derive the 
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This removes the dependence on luminosity and leaves 
only weak dependences on mass and effective tempera- 
ture, with the latter only varying by a small amount on 

the red-giant b ranch. 

As noted by IMosser et al.1 (|2010f ) . the difference be- 
tween the relation for r ed giants and the one including 
main-sequence stars by IStello et"aH (|2009f) is small but 
significant. To illustrate this difference we plot the rela- 
tions together with observations and evolutionary tracks 
in Figure |U Note that i/ max and Av for the models have 
been calculated using equations (HJ and ([3]). Compared 
to Figure we replaced the ordinate in the upper panel 
by the ratio v max / Av, which scales with the radial or- 
der of maximum power. The lower panel shows this ra- 
tio with i/ max raised to the power of 0.75, which better 
illustrates the mass dispersion since in this plot, con- 
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TABLE 1 

Coefficients of the Aj/-y max relation. 



Source 


a 





"max 


# of stars 




(//Hz) 




OHz) 




A2Z 1 


0.266±0.004 


0.759±0.003 


8-221 


810 


CAN 1 


0.286±0.004 


0.745±0.003 


12-228 


893 


COR 1 


0.279±0.004 


0.749±0.003 


2-208 


1151 


DLB 1 


0.293±0.009 


0.747±0.008 


1-192 


256 


OCTI 1 


0.273±0.004 


0.752±0.003 


14-208 


829 


OCTII 1 


0.254±0.004 


0.764±0.003 


16-203 


829 


SYD 1 


0.268±0.004 


0.758±0.003 


8-220 


799 


CoRoT 2 


0.28±0.02 


0.75±0.01 


2-100 


1827 


MS+RG 3 


0.263±0.009 


0.772±0.005 


15-4500 


55 



1 Kepler red g iants: A2Z - IMathur et aLI j20Toh. CAN - 
IKallinger eTaT] l|2010bl 'l. COR - Mosscr & Appourchaux ([200911 
DLB - Buzasi e t al. (unpublished). OCTI & OCTII - IHekker et al.l 
fl2010df) , SYD - IHuber et aLI pOOg). 

2 CoRoT red giants, see lMosser et al.l 1 12010 ). 

3 Main -sequence (MS) and red giant (RG) stars, see IStello et al.l 
(2009) and references therein. 

stant masses correspond to almost horizontal lines (see 
Equation (fTO])). The Kepler stars are shown in bins of 
20/xHz (red diamo nds). Note that the red giants used by 
IStello etldl (|2009t ) have been omitted for clarity, but do 
not differ significantly from the Kepler sample. The un- 
certainties for the main-sequence sample have been col- 
lected from the literature where available or were set to 
typical values of 3% for i/ max and 1% for Av. 

The scatter about the Av-v max relation is considerably 
larger than the measurement uncertainties (which are 
~1% for Av in our sample) and for red giants is mainly 
caused by the spread of stellar masses (Kallingc r et al.l 
l2010d) . While Figure |4] shows that this spread is signif- 
icant on the red-giant branch, it can be seen that the 
evolutionary tracks (solid lines) almost overlap for main- 
sequence stars. This effect causes a different fit to the re- 
lation between Av and v max for a sample consisting only 
of red giants compared to when main-sequence stars are 
included, which is clearly reflected by the observations 
and the determined fits shown in Figure HI 

The dotted and dashed triple-dotted lines in Figure 
[4] illustrate the effect of changing the metallicity for 
the f Mq tra ck. The effect is relat ively small on the 
giant branch (Kallingc r et al.l |2010cD . but more signifi- 
cant for less evolved stars. In fact, both main-seq uence 
stars significantly above the IStello et"aTI ( |2009t ) rela- 
tion, n Boo and /iHer, are observed t o be metal-rich 
(jCarrier et al.ll2005t lYang fc Meng|[20i0l ). while the sub- 
giant t/Indi wh i ch fal ls below the relation is metal-poor 
(|Bedding et al.l 120061) . Hence, these observations are 
qualitatively in agreement with expectations from evo- 
lutionary models combined with scaling relations. 

A detailed investigation of these effects using pulsa- 
tion models (which would be necessary to quantify mass 
and metallicity changes in Figure 3] in absolute terms) 
is beyond the scope of this paper. However, our results 
show that the observed difference in the Av-v niax rela- 
tion for red giants compared to when main-sequence stars 
are included is consistent with scaling relations. A more 
detailed investigation of the relation for main-sequence 
stars using Kepler short-cadence data will be presented 
in a forthcoming paper. 
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Fig. 4.— Upper panel: Solar-metallicity (Y=0.283, Z=0.017) 
evolutionary track s (solid lines) together with the relation by 
IStello e t al. (2009) (blue dashed-dotted line) and the relation fit- 
ted to the Kepler data (red dashed line). Blue squar e s sho w the 
main-sequence and sub-giant stars used by Stcllo ct al. (2009|), and 
red diamonds the Kepler observations as shown in Figure [3] in 
bins of 20 /tHz. The black dashed triple-dotted and dotted lines 
show the solar-mass tracks for chemical compositions of (Y,Z) = 
(0.291,0.009) and (Y,Z) = (0.265,0.035), respectively. Lower panel: 
Same as the upper panel, but with the luminosity dependence of 
the ordin ate removed by raising i/ max to the power of 0.75 (see 
Equation 1(10) '), 



4.2. Population effects in the f max and Av distributions 

The large number and diversity of red giants for which 
solar-like oscillations have been detected enables us to 
identify stellar populations using global oscillation pa- 
rameters. Due to the straightforward relationship be- 
tween oscillation parameters and fundamental parame- 
ters (see equations ([TJ) and ([3])), comparisons with mod- 
els could allow us to draw some conclusions about the 
star-formation history in the Kepler field. Such studies 
have alr eady been performed for red giants o bserved by 
CoRoT (jMiglio et al.l 12001 lYang et al1l20Tnh . 

Before such a comparison is made, it is important to 
consider potential biases in the sample of Kepler red 
giants. Firstly, limitations on the bandwidth of space- 
craft communications restrict the number of stars for 
which data can be obtained. Hence, absolute star counts 
for a given field and magnitude range will be incom- 
plete. Secondly, about two thirds of the Kepler red-giant 
sample have been chosen as astrometric reference stars 
(|Batalha et alj 120101 : iMonet et ail 120101 ). These have 
been selected to be distant, bright but unsaturated gi- 
ants (T cff ,Kic < 5400 K, loggKic < 3.8, R/R q ,kic > 2) 
in the Kepler magnitude range 11.0-12.5. Through the 
remainder of this paper we will refer to these stars as the 
astrometric sample, and to the remaining stars as the 
asteroseismic sample. Thirdly, as was the case for the 
CoRoT sample, the Kepler sample excludes stars pul- 
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Fig. 5. — Histograms of f ma x (left panels) and Au (right panels) in percent comparing a synthetic stellar population for the Kepler FOV 
(top panels), the observed distributions in the Kepler sample (middle panels) and the observed distributions in the CoRoT sample (bottom 
panels). The dashed and dotted lines in the middle panel separately show the distributions of the astrometric and asteroseismic sample, 
respectively. Note that the comparison for the Kepler sample should be restricted to stars with i/ max > lO/iHz (Ai/ > 1.5 /xHz). 



sating at frequencies that are too low to be resolved with 
the given time base available at this stage of the mission 
(^max < 10 //Hz). 

To obtain a first qualitative analysis of the observed 
distributions of the Kepler sample, we calc ulated a syn- 
thetic stellar population the same way as iMiglio et al.l 
(2009) , using the stellar synthesis code TRILEGAL 
( Girardi et al.|[2005l ). The input par ameters were i denti- 
cal to the simulations performed by Mig lio et al.l (2009) 
for the CoRoT sample, with the exception of restricting 
the simulation to a 10 deg 2 field centered on the Kepler 
FOV (a=290.7°, (5=44.5°) in the magnitude range 8-13 
mag in the Kepler bandpass. A constant star-formation 
rate was assumed. 

Figure [5] compares the distributions of the synthetic 
Kepler population w ith the observed on e, together with 
the CoRoT sample (jMosser et al.ll2010f ). Note that the 
latter includes stars from both fields of view of CoRoT, 
one directed to the galactic center, the other to the galac- 
tic anti-center. As for the CoRoT sample, we observe 
a maximum at ^ ma x^ 30 //Hz, corresponding to the red 
clump. These are low-mass (< 2Mq) stars, for which He- 
core burning occurs at similar luminosities (and hence 
similar v max ), forming a dominan t population on the 
giant branch. IMiglio et al.1 (|2009f ) noted that for the 
CoRoT stars, the maximum appears at a lower value 
of z^max (higher luminosity) than for the model. This is 
not the case for the Kepler field. 

The dashed and dotted lines in the middle panel show 
the distributions of the astrometric and asteroseismic 
sample separately, to illustrate potential biases intro- 
duced in the overall distribution by combining these dis- 
tinctly selected sets of stars. We see that the red clump 
maximum is predominantly formed by the astrometric 
sample, while the asteroseismic sample contributes more 
to slightly higher f max values. This suggests that the 



astrometric sample is largely unbiased, while the broad- 
ening of the red clump peak compared to the model is 
potentially caused by selection bias in the asteroseismic 
sample. 

Apart from the red clump, a second and much broader 
component in the synthetic population can be identified 
in the interval z^ ma x — 40—110 /-iHz and Ais ~ 5 — 10 /iHz. 
This consists of more massive stars that, compared to the 
red clump, occupy lower luminosities (and hence higher 
^max) over a wider range when they settle as He-core 
burning stars. T his component is referred to as the sec- 
ondary clump bv lGirardil (|1999| ). A comparison with the 
histogram of the Kepler sample shows that the observa- 
tions qualitatively reproduce the distributions for both 
higher v max and Az/ values. 

To investigate this further, Figure [6] shows the ratio 
z^ max /Ay as a function of i/ max for the synthetic popula- 
tion and for the Kepler observations. The solid lines are 
evolutionary tracks of different masses. Figure [7] shows 
the same plot but with the luminosity dependence on the 
ordinate removed by raising i/ max to the power of 0.75. 
In both figures, the separation of the red clump and the 
secondary clump is obvious in the population synthesis 
(upper panels). For the observations (lower panels) we 
can identify a group of stars with M > 2 M© extend- 
ing up to z^max ~ HO /iHz, as expected for the secondary 
clump population. 

The qualitative agreement between the distributions 
is encouraging and the excess of stars with i/ maj 
■l!) - 110 /iHz and M > 2 Af Q suggests that the Kepler 
sample includes stars belonging to the s e conda ry clump 
population. As pointed out by iGirardl (|1999| ). the de- 
tection of these stars should allow us to probe important 
physics, such as convective-core overshooting, and help 
to put tight constraints on the recent star-formation his- 
tory in the galaxy. The detailed modelling required for 
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Fig. 6. — Ratio of u max /Au as a function of v max for the synthetic 
population (upper panel, 1530 stars) and the Kepler red-giant 
sample (lower panel, 801 stars). Red lines show solar- metallicity 
(Y=0.283, Z=0.017) ASTEC evolutionary tracks with masses as 
indicated in the plots. The red clump and secondary clump popu- 
lations are indicated for the synthetic population. 

these inferences is beyond the scope of this paper. 

4.3. Small frequency separations 

As discussed in Section 13. 2\ the small frequency sep- 
arations of main-sequence stars depend on the sound- 
speed gradient in the stellar core and are sensitive to the 
evolutionary state of a star. The excellent quality of the 
Kepler data allows us to measure the small separations in 
an unprecedented number of stars covering a wide range 
of evolutionary states on the giant branch. 

To measure the small separations, we first examined 
the echelle diagram of each star using the power spec- 
trum in the frequency range i/ max ± 5Av and then man- 
ually fine-tuned the large separation to make the I = 
ridge vertical. The adjustment to Av was typically a few 
tenths of a microhertz. Stars for which no unambiguous 
identification of I — and 2 could be found were dis- 
carded from this analysis. The adjusted echelle diagram 
was then collapsed, and a Gaussian function was fitted 
to each mode ridge that was identified. The center of the 
fitted Gaussian was taken as the position of that ridge. 
Using this technique we were able to measure the I = 
and 2 ridges for 470 stars. Of these, we measured the 
I = 1 ridge for 400 stars and the I — 3 ridge for 45 stars. 
The uncertainties were determined from extensive simu- 
lations of artificial Kepler data that were analyzed using 
the same method. Before discussing the separations of 
the ridges, we first consider their absolute positions. 

4.3.1. Variation of e 

The measured mode ridge centroids allow us to in- 
vestigate the parameter e in Equation ([2]). Guided by 
the Sun, for which e ~ 1.5 (see, e.g., Equation (11) in 
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Fig. 7. — Same as Figure|51 but with the luminosity dependence 
of the ord inat e removed by raising f max to the power of 0.75 (see 
Equation lfTOt '1 . 



iKieldsen fc Be dding 199^), we have plotted the centroid 
positions in Figure [5J which shows that e is a function of 

Figure [8] also shows a considerable spread in e for a 
given Au, in particular for stars with Av < 5 /iHz. While 
the spread is of the same order as the estimated uncer- 
tainties, we have tested whether some of it might be re- 
lated to physical properties of the stars. We have in- 
vestigated this by correlating e with ^ max / Av which, as 
shown in Figure HI is sensitive to stellar mass. No clear 
correlation could be found. 

The relation between e and Av (and hence also 
f-max) implies that e is a function of fundamental pa- 
rameters. If this relation can be quantified for less- 
evolved stars, clear mode ridge identifications as pre- 
sented in this paper could potential ly be used to predict 
e for o ther stars. As suggested bv iBedding &; Kjeldsenl 
(2010), such comparisons can be of considerable help 
in cases where the mode ridge identification is diffi- 
cult due to short mode lifetimes or rotational splitting 
(see, e.g..lAppourchaux et al.ll200c|lBenomar et al1l2009l: 
Garcia et al.ll2009t TKallingc r et aLll2010at IBedding et al.l 
201QajT~ 

Meanwhile, we can already use the ensemble results 
in Figure |8] to suggest ridge iden tifications for the fou r 
CoRoT red giants discussed by iHekker et al.l (|2010aD . 
These stars have Av in the range 3.1 to 5.3 fiHz and so we 
expect e to be 1.0 or slightly less. Looking at the echelle 
diagrams of these st ars (see Figures 8, 10, 12 and 14 in 
IHekker et al1l2010al) . we indeed see in all cases that one 
of the two ridges falls there. In each of these diagrams, 
this indicates it is the left-hand ridge that corresponds 
to I = 1 and the other to I = 0. 

4.3.2. C-D diagrams 
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Fig. 8. — Mode ridge centroids as measured from the folded and 
collapsed power spectrum for £ = (black squares), 1 = 1 (grey 
diamonds) and I = 2 (grey triangles) . Note that the upper abscissa 
label shows e for 1 = modes. Error bars are only shown for I = 
for clarity. 

Using the measured mode ridge centroids, we were 
able to derive the small separations Svq2, Si/qi and Svq^ 
for 470, 400 and 45 stars, respec tively. The results are 
show n in so-called C-D diagrams ([Christensen-Dalsgaardl 
1988) in Figure [§1 We also plot the frequency separa- 
tion ratios Svqi/Av which, according to models, are ex- 
pected to be largely insensiti v e to surface layer effects 
([Roxburgh k. Vorontsovl l2003t lOti Floranes et all 120051; 
IMazumdarll2005l ; IChaplin et al.ll2005lT ~ 

We can see in Figure that Sv 02 is an almost 
fixed fraction of Av, which confirms the findings by 
iBedding et all (|2010bf) for low-luminosity red giants. We 
obtain the following relation for the full range of u max : 

5v 02 = (0.121 ± 0.001)Ai> + (0.047 ± 0.008) , (11) 

in agre ement with the relation found by Bedding et al. 
(2010b). We see a deviation from this linear relationship 
for stars with Av < 10/xHz (v max < 120 /iHz), which can 
be seen in Figure IHb as an increase of the ratio 61/02/ Av. 
Additionally, we observe an increased spread of 5 ^02, par- 
ticularly for stars with Av ~ 4 /xHz. While most of this 
spread in the red clump is due to the much larger number 
of stars (see Figure [5]) and the larger uncertainty in de- 
termining 5 vq2 , we tested whether it could also be related 
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Fig. 9. — C-D diagram of Su 01 versus A.v for I = 2 (a), I = 1 (c) 
and I = 3 (e) and frequency separation ratios of Svqi/Av for I = 2 
(6), I = 1 (d) and I = 3 (/). Dashed lines in panels (a) and (e) 
show linear fits to the data. 



to physical properties of the stars. C-D diagrams for red 
giants calculated from stellar models (T. R. White et al, 
in preparation) show that the expected range in &v§ij Av 
for M = 1 — 2 Mq and solar metallicity in non He-core 
burning models is about 0.02, which is roughly compara- 
ble to the range of values observed in the data outside the 
red clump. Indeed, a comparison of stars with different 
5vq2 in the lower panel of Figure [6] has shown the spread 
can be partially explained by a spread in stellar masses, 
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with lower-mass stars showing higher Svq2 values. More 
quantitative conclusions about the mass spread will be 
possible once the measurement uncertainties are reduced 
by the collection of more data. 

Figure |5t shows that 8v 01 is ne gative for almo s t all re d 
giants, confirming the findings bv: Beddi ng et al.l ((2010b). 
As for I = 2, we observe a trend (but with opposite 
sign) in the frequency separation ratio Svqi/ Av with Av, 
which can be seen in Figure[9jl. The decrease of Svoi/Av 
and the increase of bv^l Av appear to affect stars in 
the same range of Av and by a similar amount (~5%). 
Using stars common to both samples, we calculated a 
correlation coefficient of —0.1. While this indicates that 
the observed trends are statistically almost uncorrelated, 
we cannot at this point exclude that this is due to the 
large uncertain ty of measuring 5vm . 

As shown by Beddin g" et al.1 ()2010bl ) , the Kepler data 
allow the detection of I = 3 modes in red giants. We 
were able to measure the small separation Svos for 45 
stars and a linear fit to the relation with Av, shown in 
Figured, yields: 



5v m = (0.282 ± 0.005)Ai/ + (0.16 ± 0.04) 



(12) 



Figure [9f shows the ratio Svq^/ Av, which shows some ev- 
idence for an increase of Svqs with decreasing Av. While 
the error bars appear to be too large to make any firm 
statements on this variation, the observed scatter sug- 
gests that the uncertainties are overestimated. Given 
that the relative amplitude of Z = 3 modes is not well 
known, it seems likely that the strength of I = 3 modes 
has been considerably underestimated in the simulations 
used to derive the uncertainties. 

4.3.3. Ensemble collapsed echelle diagram 

In order to display the separations of the mode ridges 
for all stars, we use the measured values of e to shift 
and align the folded and collapsed power spectra in an 
ensemble collapsed echelle diagram. The result is shown 
in the upper panel of Figure 1101 Note that this diagram 
is slightly different from the sca led echelle diagram pre- 
sented by Bedding et al. (2010b) in that we have shifted 
the folded and collapsed power spectra rather than scal- 
ing frequencies. The thick solid line in the lower panel 
shows the upper panel summed along the full range 
of z^max- It clearly reveals the presence of ridges with 
/ = 0, 1, 2, and also 1 = 3. Th e power for I = 3 is weaker 
than found by Bedding et al. (2010b), who considered a 
smaller sample of low-luminosity giants with high signal- 
to-noise. We attribute this to the larger (and less biased) 
sample of stars in this paper, including stars with unde- 
tectable I = 3 modes, as well as to the fact that the I = 3 
ridge has a slight tilt, making its power spread out in the 
lower panel of Figure [TO] 

The trends of the frequency separation ratios observed 
in Figure [9] are clearly visible in Figure [101 As i/ max (and 
hence Av) decreases, the I = 2 ridge slopes to the left 
while the I = 1 ridge slopes to the right, corresponding to 
an increase of S v^l Av and a decrease of Sv^i/ Av. These 
features are also clearly shown by the thin red and blue 
lines in the lower panel of Figure [TUJ which correspond 
to subsets with f max < 50/iHz and v max > 100 /iHz, 
respectively. The red and blue lines also show evidence 
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Fig. 10. — Upper panel: Ensemble collapsed echelle diagram of 
470 red giants. Note that the contour threshold has been lowered to 
40% of the maximum value and that each row shows the folded and 
collapsed power spectrum of one star. The right ordinate marks 
Vmax values for selected rows. Lower panel: Upper panel collapsed 
along the entire range of v maI (thick black line) as well as in subsets 
i/ msI > 100 /^Hz (thin blue line) and ^ max < 50 /iHz (thin red line). 
Ridge identification and definitions of small separations used in 
this paper are indicated. In both panels, the dotted lines mark the 
centers of the I = 0, 2 and 3 ridge and the dashed line shows the 
midpoint of adjacent I = modes. 



for a slope of the / = 3 ridge, corresponding to an increase 
of 5vo3 as observed in Figure [9f. 

We note that the relative width of the I = ridge 
shown in Figure fTOl (measured in terms of A^) decreases 
with increasing f max . In fact, the absolute width of the 
ridge (in /iHz), which is determined by the frequency res- 
olution, mode lifetime and curvature, remains roughly 
co nstant over the range of stars considered. As found 
bv lBedding et all (|2010b| ) for low-luminosity red giants, 
the I = 1 ridge is significantly broader than the others. 
This was interpr eted as a dire c t con firmation of theo- 
retical results by iDupret et al.l (|2009f) , which predicted 
complicated power spectra due to less efficient trapping 
of mixed modes in the cores of low-luminosity red giants. 
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Fig. 11. — Absolute width of the I = (dashed line), I = 1 (solid 
line) and I = 2 (dashed-dotted line) ridge as a function of v max . 
The dotted line marks the formal frequency resolution of the data 
(0.08 fiUz). 

The models by iDupret et al.l ()2009[ ) also predicted that 
this effect should become less pronounced for higher lu- 
minosity red giants. Indeed, we observe that the relative 
width of the I = 1 ridge shown in Figure [TU] remains 
roughly constant, indicating that its absolute width is 
significantly lower for high-luminosity stars than for low- 
luminosity stars. 

To quantify this, we measured the absolute widths of 
the / = 0, 1 and 2 ridges in several bins of f ma x, and show 
the result in Figure [TU Note that we have accounted for 
possible artificial broadening of the ridges due to Nyquist 
effects by excluding the two stars in our sample with 
i/ max > 230/iHz from this calculation. As can be seen in 
Figure [TTJ the widths of the I = and 2 ridges remain 
roughly constant over the range of ^ max , while the width 
of the I = 1 ridge increases significantly. Our observation 
of a much narrower I = 1 ridge for high-luminosity stars 
therefore provides f urther confirmation of the results by 
IDupret et all (|2009T ). 

4.4. Pulsation amplitudes 

Amplitudes of solar-like oscillations provide valuable 
information to test models of convection (see Section 
13.31) . Equations ([TJ , ([7]) and ([8]) , in combination with ef- 
fective temperatures, can be used to t est theoretical sca l- 
ing relations for amplitudes (see, e.g-. IStello et al.| [2007). 
as follows: 

As shown bv lMosser et al.l (|2010l ) for the CoRoT red gi- 
ants, the relation between f max and the pulsation ampli- 
tudes indeed follows a power law and this is shown for 
our Kepler sample in Figure 1121 We have investigated 
the slope of the power law using amplitudes from four 
different methods (A2Z, COR, OCT and SYD, see Table 
1 for references) and find on average: 

A 650n m oc v m °J . (14) 

As can be seen from Figure [T2] the measured ampli- 
tudes show considerable structure, in particular a lack 
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Fig. 12. — Pulsation amplitude versus f ma x observed in the 
Kepler bandpass (A = 650 nm). The red solid line shows a power 
law using a slope of —0.8. Red symbols show all stars with 
fmax > 40 /iHz and M > 2Mq identified in Figures \E\ and [7] 



of low amplitudes for i/ max > 110 /iHz. This feature is 
found consistently in all methods and therefore appears 
to be intrinsic. Some low-amplitude stars do appear at 
f-max si 180 /iHz, but we note that these amplitude esti- 
mates are rather uncertain due to the difficulty of esti- 
mating the background noise for stars oscillating near the 
Nyquist frequency. To investigate this lack of low ampli- 
tudes, we plot all stars that we have tentatively identified 
as secondary clump stars (using Figures [5] and setting 
f-Wx > 40 /iHz and M > 2M Q ) in Figure [12] as red sym- 
bols. We observe that these stars have systematically low 
amplitudes. We conclude that, for a given ^ max in the 
range between 40 — 110 /iHz, stars with lower Av (and 
therefore higher mass, see Equation (10)) oscillate with 
lower amplitudes than other stars with the same i/ max . 
Assuming that this also holds for low-luminosity red gi- 
ants, the lack of low amplitude stars for ^ max > 110 /iHz 
could then possibly be explained by the relatively fast 
evolution of more massive stars near the bottom of the 
red-giant branch compared to ~ 1 Mq stars. 

Both the observation of low amplitudes for stars with 
lower A.v in the range f max — 40 — 110 //Hz and the lack 
of low amplitude stars for f max > 110 /iHz indicate a 
significant dependence of the amplitude- ^ max relation on 
the stellar mass. Since such a dependence is not seen 
in Equation (fT3"]l . this suggests a revision of the scaling 
relation in Equation ([7]). 

The overall distribution of the measured amplitudes in 
Figure [12] shows considerable additional scatter, which 
is possibly related to p hysical effects such as metallic- 
ity (|Samadi et al. 2010a;IJ) in combination with measure- 
ment uncertainties. Tests have shown that the measure- 
ment error for amplitudes is dominated by the specific 
method and model used to subtract the background sig- 
nal due to granulation and stellar activity. Further work 
is needed to identify possible systematic effects for the 
entire range of ^ max , which will then allow us to improve 
the relation between the oscillation amplitude and ^ max 
using a large sample of stars on the giant branch. 

5. SUMMARY AND CONCLUSIONS 
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Using data from the first four months of the Kepler 
mission we have studied various aspects of solar-like os- 
cillations in ^800 red giants. Our conclusions can be 
summarized as follows: 

1. The Av-v mlix relation using red giants is signif- 
icantly different from the relation when main- 
sequence stars are included. This difference can 
be explained by evolutionary models and scaling 
relations, which show that the Ai^-^ max relation is 
mainly sensitive to stellar masses for red giants and 
to metallicity for main-sequence stars. 

2. The v max and Av distributions are in qualitative 
agreement with a simple model of the stellar pop- 
ulation. We observed that v max for the red clump 
stars in the Kepler field is slightly higher than for 
the CoRoT fields, possibly due to the different 
stellar populations that are observed. We identi- 
fied several stars with v max = 40 — 110 /iHz and 
M > 2M Q , in agreement with a secondary clump 
population. 



evidence that this is partially due to the spread in 
mass. We also measured the small separation 8v^ 
and tentatively identify a similar variation of the 
frequency separation ratio 5 v 03 / 'Av as observed for 
5vo2/ 'Av as a function of Av. 

5. The absolute width of the I = 1 ridge for stars of 
higher luminosity is significantly narrower than for 
stars with low luminosity. This is the first quanti- 
tative confirmation of more efficient mode trapping 
as predicted by theory. 

6. We presented a first estimate for the relation be- 
tween pulsation amplitude and v max for Kepler 
red giants. We observed a distinct lack of low- 
amplitude stars for ^ max > 110 /iHz and found that, 
for a given v max , stars with lower Av (and therefore 
higher mass, see Equation (10)) tend to show lower 
amplitudes than stars with higher Av. Both obser- 
vations can be explained with a mass dependence of 
the amplitude-j/ max relation, and therefore suggest 
that the scaling relation for luminosity amplitudes 
of red giants needs to be revised. 



3. The quantity e in the asymptotic relation is a func- 
tion of Av. If the relation between e and fundamen- 
tal parameters can be quantified and is also valid 
for main-sequence stars, it could be used to study 
surface layer effects or for mode identification in 
stars where ridges cannot be clearly identified due 
to short mode lifetimes or rotational splitting. We 
demonstrated the potential of the observed rela- 
tion by providing mode identifications for four red 
giants observed by C0R0T. 

4. We presented C-D diagrams for I = 1, 2 and 3 and 
show that frequency separation ratios of S V02/ Av 
and Svqi/Av reveal opposite trends as a function of 
Av. We observed a spread in the small separations, 
in particular for stars in the red clump, and find 
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